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Abstract 

We introduce a set of correlations between energy flow and event shapes that are sensitive 
to the flow of color at short distances in jet events. These correlations are formulated for 
a general set of event shapes, which includes jet broadening and thrust as special cases. 
We illustrate the method for e + e~ dijet events, and calculate the correlation at leading 
logarithm in the energy flow and at next-to-leading-logarithm in the event shape. 



1 Introduction 

The agreement of theoretical predictions with experiment for jet cross sections is often impres- 
sive. This is especially so for inclusive jet cross sections at high px, using fixed-order factorized 
perturbation theory and parton distribution functions [1]. A good deal is also known about the 
substructure of jets, through the theoretical and experimental study of multiplicity distributions 
and fragmentation functions [2], and of event shapes [3, 4, 5]. Event shape distributions [6, 7, 8] 
in particular offer a bridge between the perturbative, short-distance and the nonperturbative, 
long-distance dynamics of QCD [9]. 

Energy flow [10] into angular regions between energetic jets gives information that is in some 
ways complementary to what we learn from event shapes. In perturbation theory, the distribution 



1 



of particles in the final state reflects interference between radiation from different jets [2], and 
there is ample evidence for perturbative antenna patterns in interjet radiation at both e + e~ [11] 
and hadron colliders [12, 13]. Energy flow between jets must also encode the mechanisms that 
neutralize color in the hadronization process, and the transition of QCD from weak to strong 
coupling. Knowledge of the interplay between energy and color flows [14, 15] may help identify 
the underlying event in hadron collisions [16], to distinguish QCD bremsstrahlung from signals 
of new physics. Nevertheless, the systematic computation of energy flow into interjet regions 
has turned out to be subtle [17] for reasons that we will review below, and requires a careful 
construction of the class of jet events. It is the purpose of this work to provide such a construction, 
using event shapes tool. 

In this paper, we introduce correlations between event shapes and energy flow, "shape/flow 
correlations" , that are sensitive primarily to radiation from the highest-energy jets. So long as the 
observed energy is not too small, in a manner to be quantified below, we may control logarithms 
of the ratio of energy flow to jet energy [15, 18]. 

The energy flow observables that we discuss below are distributions associated with radiation 
into a chosen interjet angular region, Q. Within fl we identify a kinematic quantity Qq = eQ, 
at cm. energy Q, with £ < 1. Qq may be the sum of energies, transverse energies or related 
observables for the particles emitted into fi. Let us denote by the complement of f2. We are 
interested in the distribution of Qq, for events with a fixed number of jets in Cl. This set of events 
may be represented schematically as 

A + B -> Jets + X n + R n (Qn) ■ (1) 

Here X n stands for radiation into the regions between Q and the jet axes, and Rq for radiation 
into Q. 

The subtlety associated with the computation of energy flow concerns the origin of logarithms, 
and is illustrated by Fig. 1. Gluon 1 in Fig. 1 is an example of a primary gluon, emitted directly 
from the hard partons near a jet axes. Phase space integrals for primary emissions contribute 
single logarithms per loop: (I/Qq)^ \n n ~ 1 (Q/Q n ) = (l/eQ)a^\ii n ^ l (l/e), n > 1, and these 
logarithms exponentiate in a straightforward fashion [15]. At fixed Qq for Eq. (1), however, 
there is another source of potentially large logarithmic corrections in Q n . These are illustrated 
by gluon 2 in the figure, an example of secondary radiation in Q, originating a parton emitted by 
one of the leading jets that define the event into intermediate region Cl. As observed by Dasgupta 
and Salam [17], emissions into Q from such secondary partons can also result in logarithmic 
corrections, of the form (l/Qnja™ ln n ~ 1 (Qn/Qn), n > 2, where is the maximum energy 
emitted into Cl. These logarithms arise from strong ordering in the energies of the primary and 
secondary radiation because real and virtual enhancements associated with secondary emissions 
do not cancel each other fully at fixed Qq. 

If the cross section is fully inclusive outside of Q, so that no restriction is placed on the 
radiation into 0, Qq can approach Q, and the secondary logarithms can become as important 
as the primary logarithms. Such a cross section, in which only radiation into a fixed portion of 
phase space is specified, was termed "non-global" by Dasgupta and Salam, and the associated 
logarithms are also called non-global [17, 19, 20]. 
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In effect, a non-global definition of energy flow is not restrictive enough to limit final states 
to a specific set of jets, and non-global logarithms are produced by jets of intermediate energy, 
emitted in directions between region Q and the leading jets. Thus, interjet energy flow does 
not always originate directly from the leading jets, in the absence of a systematic criterion for 
suppressing intermediate radiation. Correspondingly, non-global logarithms reflect color flow at 
all scales, and do not exponentiate in a simple manner. Our aim in this paper is to formulate a 
set of observables for interjet radiation in which non-global logarithms are replaced by calculable 
corrections, and which reflect the flow of color at short distances. By restricting the sizes of event 
shapes, we will limit radiation in region Q, while retaining the chosen jet structure. 

An important observation that we will employ below is that non-global logarithms are not 
produced by secondary emissions that are very close to a jet direction, because a jet of parallel- 
moving particles emits soft radiation coherently. By fixing the value of an event shape near 
the limit of narrow jets, we avoid final states with large energies in Q away from the jet axes. 
At the same time, we will identify limits in which non-global logarithms reemerge as leading 
corrections, and where the methods introduced to study nongobal effects in Refs. [17, 19, 20] 
provide important insights. 




Figure 1: Sources of global and non-global logarithms in dijet events. Configuration 1, a primary 
emission, is the source of global logarithms. Configuration 2 can give non-global logarithms. 

To formalize these observations, we study below correlated observables for e + e~ annihilation 
into two jets. (In Eq. (1) A and B denote positron and electron.) In e + e~ annihilation dijet 
events, the underlying color flow pattern is simple, which enables us to concentrate on the energy 
flow within the event. We will introduce a class of event shapes, f(a) suitable for measuring 
energy flow into only part of phase space, with a an adjustable parameter. To avoid large 
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non-global logarithmic corrections we weight events by exp[— uf], with v the Laplace transform 
conjugate variable. 

For the restricted set of events with narrow jets, energy flow is proportional to the lowest-order 
cross section for gluon radiation into the selected region. The resummed cross section, however, 
remains sensitive to color flow at short distances through anomalous dimensions associated with 
coherent interjet soft emission. In a sense, our results show that an appropriate selection of jet 
events automatically suppresses nonglobal logarithms, and confirms the observation of coherence 
in interjet radiation [2, 12]. 

In the next section, we introduce the event shapes that we will correlate with energy flow, 
and describe their relation to the thrust and jet broadening. Section 3 contains the details of the 
factorization procedure that characterizes the cross section in the two-jet limit. This is followed 
in Sec. 4 by a derivation of the resummation of logarithms of the event shape and energy flow, 
following the method introduced by Collins and Soper [21]. We then go on in Sec. 5 to exhibit 
analytic results at leading logarithmic accuracy in Qn/Q and next-to- leading logarithm in the 
event shape. Section 6 contains representative numerical results. We conclude with a summary 
and a brief outlook on further applications. 

2 Shape/Flow Correlations 

2.1 Weights and energy flow in dijet events 

In the notation of Eq. (1), we will study an event shape distribution for the process 

e + + e~ - JxfrjJ + J 2 (p j2 ) + X n (/) + R n (Q n ) , (2) 

at cm. energy Q 3> Qq 3> A QCD . Two jets with momenta pj c , c = 1, 2 emit soft radiation (only) 
at wide angles. Again, Q is a region between the jets to be specified below, where the total energy 
or the transverse energy Qq of the soft radiation is measured, and Cl denotes the remaining phase 
space (see Fig. 1). Radiation into f2 is constrained by event shape /. We refer to cross sections 
at fixed values (or transforms) of / and Qq as shape/flow correlations. 

To impose the two-jet condition on the states of Eq. (2) we choose weights that suppress 
states with substantial radiation into away from the jet axes. We now introduce a class of 
event shapes /, related to the thrust, that enforce the two-jet condition in a natural way. 

These event shapes interpolate between and extend the familiar thrust [4] and jet broadening 
[7, 8], through an adjustable parameter a. For each state A" that defines process (2), we separate 
f2 into two regions, C , c = 1,2, containing jet axes, n c (N). To be specific, we let Cli and 
2 be two hemispheres that cover the entire space except for their intersections with region Q. 
Region Cli is centered on hi, and fi 2 is the opposite hemisphere. We will specify the method 
that determines the jet axes hi and n 2 momentarily. To identify a meaningful jet, of course, the 
total energy within Cli should be a large fraction of the available energy, of the order of Q/2 in 
dijet events. In e + e~ annihilation, if there is a well-collimated jet in Cli with nearly half the total 
energy, there will automatically be one in fl 2 . 
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We are now ready to define the contribution from particles in region fl c to the a-dependent 
event shape, 

fn c (N,a) = ±= £ kl ± u}- a {l-h i .h c f- a , (3) 

where a is any real number less than two, and where y/s = Q is the cm. energy. The sum is over 
those particles of state N with direction hi that flow into C , and their transverse momenta ki t ± 
are measured relative to h c . The jet axis hi for jet 1 is identified as that axis that minimizes the 
specific thrust-related quantity /^(A, a = 0). When Cl c in Eq. (3) is extended to all of phase 
space, the case a = is then essentially 1 — T, with T the thrust, while a = 1 is related to the 
jet broadening. 

Any choice a < 2 in (3) specifies an infrared safe event shape variable, because the contribution 
of any particle i to the event shape behaves as Q\~ a in the collinear limit, 9i = cos _1 (nj -h c ) — > 0. 
Negative values of a are clearly permissible, and the limit a — > — oo corresponds to the total cross 
section. At the other limit, the factorization and resummation techniques that we discuss below 
will apply only to a < 1. For a > 1, contributions to the event shape (3) from energetic particles 
near the jet axis are generically larger than contributions from soft, wide-angle radiation, or equal 
for a = 1. When this is the case, the analysis that we present below must be modified, at least 
beyond the level of leading logarithm [8] . 

In summary, once h\ is fixed, we have divided the phase space into three regions: 

• Region Q, in which we measure, for example, the energy flow, 

• Region fli, the entire hemisphere centered on h±, that is, around jet 1, except its intersection 
with Q, 

• Region 2 , the complementary hemisphere, except its intersection with fl. 
In these terms, we define the complete event shape variable / (N, a) by 

f(N,a) = f ni (N : a) + fn 2 (N,a), (4) 

with fn c , c = 1, 2 given by (3) in terms of the axes hi of jet 1 and h 2 of jet 2. We will study the 
correlations of this set of event shapes with the energy flow into f2, denoted as 

/(AO = 4 E"i- (5) 
V 6 f H en 

The differential cross section for such dijet events at fixed values of / and / is now 

wf 1 = ^ElMWiwv-P,) 

x6(e - /(JV)) 6(1 - f(N,a)) } 2 (r h - ft(JV)) , (6) 

where we sum over all final states that contribute to the weighted event, and where M(N) 
denotes the corresponding amplitude for e + e~ — > N. The total momentum is pi, with p] = s = 
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Q 2 . As mentioned in the introduction, for much of our analysis, we will work with the Laplace 
transform of (6), 



da(e,v,s,a) f°° _ vg da(e,e,s,a) 
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Jo 



de dh\ Jo de de dh\ 
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2s 



= h E \M(N)\ 2 e-^ N ^ (2n) 4 5 4 ( Pl -p N ) 

N 



x5(e-f(N))5 2 (n 1 -n(N)). (7) 

Singularities of the form (1/e) ln n (l/5) in the cross section (6) give rise to logarithms ln n+1 v in 
the transform (7). 

Since we are investigating energy flow in two-jet cross sections, we fix the constants e and e 
to be both much less than unity: 

0<£,£<1. (8) 

We refer to this as the elastic limit for the two jets. In the elastic limit, the dependence of 
the directions of the jet axes on soft radiation is weak. We will return to this dependence 
below. Independent of soft radiation, we can always choose our coordinate system such that the 
transverse momentum of jet 1 is zero, 

PJi,± = 0, (9) 

with p Jl in the x 3 direction. In the limit e, e — > 0, and in the overall cm., p Jl and p j2 then 
approach light-like vectors in the plus and minus directions: 

pj, - (vf'°"' 0j 

v% - (o + >v!' o± )- ( io ) 

As usual, it is convenient to work in light-cone coordinates, p M = (p + ,p~ ,p±), which we normalize 
as p^ = (l/y/2)(p° ±p 3 ). For small e and e, the cross section (6) has corrections in ln(l/e) and 
ln(l/e), which we will organize in the following. 



2.2 Weight functions and jet shapes 

In Eq. (3), a is a parameter that allows us to study various event shapes within the same 
formalism; it helps to control the approach to the two-jet limit. As noted above, a < 2 for 
infrared safety, although the factorization that we will discuss below applies beyond leading 
logarithm only to 1 > a > —oo. A similar weight function with a non-integer power has been 
discussed in a related context for 2 > a > 1 in [22]. To see how the parameter a affects the shape 
of the jets, let us reexpress the weight function for jet 1 as 

fd 1 (N,a) = ^= Yl ^sin a i (l-cos0 i ) 1 " a , (11) 
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where 9i is the angle of the momentum of final state particle % with respect to jet axis h\. As 
a — > 2 the weight vanishes only very slowly for 6>j — > 0, and at fixed /q i; the jet becomes very 
narrow. On the other hand, as a — > — oo, the event shape vanishes more and more rapidly in 
the forward direction, and the cross section at fixed becomes more and more inclusive in the 
radiation into 

In this paper, as in Ref. [15], we seek to control corrections in the single-logarithmic variable 
a s (Q) ln(l/e), with e = Qn/Q. Such a resummation is most relevant when 

*W)l.())>l-<«p(^;). (12) 

Let us compare these logarithms to non-global effects in shape/flow correlations. At v — and 
for a — > — oo, the cross section becomes inclusive outside Q. As we show below, the non-global 
logarithms discussed in Refs. [15, 17] appear in shape/flow correlations as logarithms of the form 
a s (Q) \n(l/(eu)), with v the moment variable conjugate to the event shape. To treat these 
logarithms as subleading for small e and (relatively) large v, we require that 

asW)ln (_L) <1 ^ E> I exp (_±_). (13) 

For large v, there is a substantial range of e in which both (12) and (13) can hold. When v is 
large, moments of the correlation are dominated precisely by events with strongly two-jet energy 
flows, which is the natural set of events in which to study the influence of color flow on interjet 
radiation. (The peak of the thrust cross section is at (1 — T) of order one-tenth at LEP energies, 
corresponding to v of order ten, so the requirement of large v is not overly restrictive.) In the next 
subsection, we show how the logarithms of (eu)^ 1 emerge in a low order example. This analysis 
also assumes that a is not large in absolute value. The event shape at fixed angle decreases 
exponentially with a, and we shall see that higher-order corrections can be proportional to a. 
We always treat In v as much larger than \a\. 

2.3 Low order example 

In this section, we check the general ideas developed above with the concrete example of a two- 
loop cross section for the process (2). This is the lowest order in which a non-global logarithm 
occurs, as observed in [17]. We normalize this cross section to the Born cross section for inclusive 
dijet production. A similar analysis for the same geometry has been carried out in [17] and [23]. 

— * 

The kinematic configuration we consider is shown in Fig. 2. Two fast partons, of velocities f3± 
and /?2, are treated in eikonal approximation. In addition, gluons are emitted into the final state. 
A soft gluon with momentum k is radiated into region Q and an energetic gluon with momentum 
I is emitted into the region 0. We consider the cross section at fixed energy, oo k = ey/s. As 
indicated above, non-global logarithms arise from strong ordering of the energies of the gluons, 
which we choose as io\ ^> uik- In this region, the gluon / plays the role of a "primary" emission, 
while A; is a "secondary" emission. 
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Figure 2: A kinematic configuration that gives rise to the non-global logarithms. A soft gluon 
with momentum k is radiated into the region Q, and an energetic gluon with momentum I is 
radiated into fl. Four- vectors (5\ and fa, define the directions of jet 1 and jet 2, respectively. 



For our calculation, we take the angular region Vt to be a "slice" or "ring" in polar angle of 
width 25, or equivalently, (pseudo) rapidity interval (—77,77), with 

A, = 2, = ln(i±^, (14) 
y 1 — sin J 

The lowest-order diagrams for this process are those shown in Fig. 3, including distinguishable 
diagrams in which the momenta k and I are interchanged. 

The diagrams of Fig. 3 give rise to color structures C% and CfCa, but terms proportional to 
Cp may be associated with a factorized contribution to the cross section, in which the gluon k is 
emitted coherently by the combinations of the gluon I and the eikonals. To generate the CfCa 
part, on the other hand, gluon k must "resolve" gluon I from the eikonal lines, giving a result 
that depends on the angles between I and the eikonal directions. 

The computation of the diagrams is outlined in Appendix A; here we quote the results. We 
adopt the notation q = cos 81, sj = sin^, with Q\ the angle of momentum I measured relative 
to fa, and similarly for k. We take, as indicated above, a Laplace transform with respect to 
the shape variable, and identify the logarithm in the conjugate variable v. We find that the 
logarithmic CVCyi-dependence of Fig. 3 may be written as a dimensionless eikonal cross section 
in terms of one energy and two polar angular integrals as 

1 \ 1 



CVCU p - dc k / dc, m-e-» u .d-^-f/O 

de \ 7T / e J-sinS JsinS Je^/s 10\ 



X 



Cfc + ql + CfcU + q 1-Cfe/ sil + ci 



(15) 



In this form, the absence of collinear singularities in the CfCa term at cos Q\ = +1 is manifest, 
independent of v. Collinear singularities in the I integral completely factorize from the k integral, 
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a) b) c) 




d) e) f) 

Figure 3: The relevant two-loop cut diagrams corresponding to the emission of two real gluons 
in the final state contributing to the eikonal cross section. The dashed line represents the final 
state, with contributions to the amplitude to the left, and to the complex conjugate amplitude 
to the right. 
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and are proportional to C\. The logarithmic dependence on e for v > 1 is readily found to be 



where C(Ar)) is a finite function of the angle S, given explicitly in Appendix A. 

We can contrast this result to what happens when v — 0, that is, for an inclusive, non-global 
cross section. In this case, recalling that e = Qci/Q, we find in place of Eq. (16) the non-global 
logarithm 



As anticipated, the effect of the transform is to replace the non-global logarithm in Q/Qn, by a 
logarithm of l/{ev). We are now ready to generalize this result, starting from the factorization 
properties of the cross section near the two-jet limit. 

3 Factorization of the Cross Section 

In this section we study the factorization of the correlations (6). The analysis is based on a 
general approach that begins with the all-orders treatment of singularities in perturbative cross 
sections [24, 25], and derives factorization from the analyticity and gauge properties of high 
energy Green functions and cross sections [26]. The functions that appear in factorized cross 
sections are expressible in terms of QCD matrix elements [27], and the matrix elements that we 
will encounter are familiar from related analyses for heavy quark and jet production [28]. We 
refer in several places below to standard arguments discussed in more detail in [25, 26]. The 
aim of this section, and the reason why a careful analysis is necessary, is to identify the specific 
dimensionless combinations of kinematic variables on which the factorized matrix elements may 
depend. We will use these dependences in the following section, when we discuss the resummation 
properties of our correlations. 

3.1 Leading regions near the two-jet limit 

In order to resum logarithms of e and e (or equivalently u, the Laplace conjugate of e) we have 
first to identify their origin in momentum space when e, e — > 0. Following the procedure and 
terminology of [24], we identify "leading regions" in the momentum integrals of cut diagrams, 
which can give rise to logarithmic enhancements of the cross section associated with lines ap- 
proaching the mass shell. Within these regions, the lines of a cut diagram fall into the following 
sub diagrams: 

• A hard-scattering, or "short-distance" subdiagram H, where all components of line mo- 
menta are far off-shell, by order Q. 




(16) 




(17) 
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• Jet sub diagrams, J\ and J2, where energies are fixed and momenta are collinear to the 
outgoing primary partons and the jet directions that emerge from the hard scattering. (For 
e = e = 0, the sum of all energies in each jet is one-half the total energy.) To characterize 
the momenta of the lines within the jets, we introduce a scaling variable, A < 1. Within 
jet 1, momenta i scale as (i + ~ Q,i~ ~ \Q,(± ~ A 1 / 2 ^). 

• A soft subdiagram, S connecting the jet functions J\ and J 2 , in which the components of 
momenta k are small compared to Q in all components, scaling as (k^ 1 ~ \Q, k± ~ XQ). 

An arbitrary final state N is the union of substates associated with these subdiagrams: 

N = N s © Nj, © N J2 . (18) 

As a result, the event shape / can also be written as a sum of contributions from the soft and 
jet subdiagrams: 

f(N,a) = f N (N s ,a) + fg(N h ,a) + fl(N h ,a). (19) 

The superscript N reminds us that the contributions of final-state particles associated with the 
soft and jet functions depend implicitly on the full final state, through the determination of the 
jet axes, as discussed in Sec. 2. In contrast, the energy flow weight, f(N), depends only on 
particles emitted at wide angles, and is hence insensitive to collinear radiation: 

f(N) = f(N s ) . (20) 

When we sum over all diagrams that have a fixed final state, the contributions from these 
leading regions may be factorized into a set of functions, each of which corresponds to one of the 
generic hard, soft and jet subdiagrams. The arguments for this factorization at leading power 
have been discussed extensively [21, 26, 29]. The cross section becomes a convolution in e, with 
the sums over states linked by the delta function which fixes hi, and by momentum conservation, 



da(e, e, s, a) dcr 
dededhi dh 



^H{sM) E [de s S(N s )5(e-f(N s ))6(e s -f N (N s ,a)) 

Ns,N Jc J 

f[ / de Jc J c {N Jc ) 5{e Jc - fl(N Jc , a)) 

c=l J 



x (2^5^p I -p(Nj 2 )-p(N Jl )-p(N s )) 
x 5 2 (h 1 - n(N)) 5{e-e, h - e, h - e s ) 

dao 



dh\ 



5(e) 5(e) +<D(a s ). (21) 



Here da^/dhi is the Born cross section for the production of a single particle (quark or anti- 
quark) in direction hi, while the short-distance function H(s,h\) = 1 + G(a s ), which describes 
corrections to the hard scattering, is an expansion in a s with finite coefficients. The functions 
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j7~ c (iVj c ), S(N S ) describe the internal dynamics of the jets and wide-angle soft radiation, re- 
spectively. We will specify these functions below. We have suppressed their dependence on a 
factorization scale. Radiation at wide angles from the jets will be well-described by our soft 
functions S(N S ), while we will construct the jet functions J c (Nj c ) to be independent of e, as in 
Eq. (21). 

So far, we have specified our sums over states in Eq. (21) only when all lines in N s are 
soft, and all lines in N Jc have momenta that are collinear, or nearly collinear to pj c . As e and 
e vanish, these are the only final-state momenta that are kinematically possible. Were we to 
restrict ourselves to these configurations only, however, it would not be straightforward to make 
the individual sums over N s and Nj c infrared safe. Thus, it is necessary to include soft partons 
in N s that are emitted near the jet directions, and soft partons in the N Jc at wide angles. We 
will show below how to define the functions J c (N Jc ), S(N S ) so that they generate factoring, 
infrared safe functions that avoid double counting. We know on the basis of the arguments of 
Refs. [21, 26, 29] that corrections to the factorization of soft from jet functions are suppressed 
by powers of the weight functions e and/or e. 

3.2 The factorization in convolution form 

Although formally factorized, the jet and soft functions in Eq. (21) are still linked in a potentially 
complicated way through their dependence on the jet axes. Our strategy is to simplify this 
complex dependence to a simple convolution in contributions to e, accurate to leading power in 
e and e. 

First, we note that the cross section of Eq. (21) is singular for vanishing e and e, but is a 
smooth function of s and hi. We may therefore make any approximation that changes s and/or 
hi by an amount that vanishes as a power of e and e in the leading regions. 

Correspondingly, the amplitudes for jet c are singular in ej c , but depend smoothly on the jet 
energy and direction, while the soft function is singular in both e and e s , but depends smoothly 
on the jet directions. As a result, at fixed values of e and e we may approximate the jet directions 
and energies by their values at e = e = in the soft and jet functions. 

Finally, we may make any approximation that affects the value of e and/or ej c by amounts 
that vanish faster than linearly for e — > 0. It is at this stage that we will require that a < 1. 

With these observations in mind, we enumerate the replacements and approximations by 
which we reduce Eq. (21), while retaining leading-power accuracy. 

1. To simplify the definitions of the jets in Eq. (21), we make the replacements f£ (Nj c , a) — > 
f c (N Jc ,a) with 

fc(N Jc ,a) = ^= £ kl ± ul- a (l-h i -h c ) 1 - a . (22) 

V s all hieN Jc 

The jet weight function f c (N Jc , a) now depends only on particles associated with N Jc . The 
contribution to f c (Nj c , a) from particles within region Cl c , is exactly the same here as in the 
weight (3), but we now include particles in all other directions. In this way, the independent 
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sums over final states of the jet amplitudes will be naturally infrared safe. The value of 
f c (Nj c ,a) differs from the value of f$ (Nj c ,a), however, due to radiation outside Cl c , as 
indicated by the new subscript. This radiation is hence at wide angles to the jet axis. In 
the elastic limit (8), it is also constrained to be soft. Double counting in contributions 
to the total event shape, f(N,a), will be avoided by an appropriate definition of the soft 
function below. The sums over states are still not yet fully independent, however, because 
the jet directions h c still depend on the full final state N. 

2. Next, we turn our attention to the condition that fixes the jet direction hi. Up to corrections 
in the orientation of hi that vanish as powers of e and e, we may neglect the dependence 
of hi on N s and Nj 2 : 

Sihi-niN^^Sihi-niNj,)). (23) 

In Appendix B, we show that this replacement also leaves the value of e unchanged, up 
to corrections that vanish as e 2 ~ a . Thus, for a < 1, (23) is acceptable to leading power. 
For a < 1, we can therefore identify the direction of jet 1 with hi. These approximations 
simplify Eq. (21) by eliminating the implicit dependence of the jet and soft weights on the 
full final state. We may now treat hi as an independent vector. 

3. In the leading regions, particles that make up each final-state jet are associated with states 
Nj c , while N s consists of soft particles only. In the momentum conservation delta function, 
we can neglect the four-momenta of lines in N s , whose energies all vanish as e, e — > 0: 

6 4 ( P i -p(N j2 ) -piNjJ -p(N s )) - 5\ PI -p j2 - PJl ). (24) 

4. Because the cross section is a smooth function of the jet energies and directions, we may 
also neglect the masses of the jets within the momentum conservation delta function, as in 
Eq. (10). In this approximation, we derive in the cm., 

^{Pi -PJ2-PJ1) <J("\/s - ^(iVjJ - u(Nj 2 )) SQpj^ - \pj 2 \) -r^-p 5 2 (hi + h 2 ) 

- -J^- "Wi)) S " ^)) ^ + fc) • (25) 
Our jets are now back-to-back: 

h 2 — ► —hi . (26) 

Implementing these replacements and approximations for a < 1, we rewrite the cross section 
Eq. (21) as 



da(e,e,s,a) da r - _ 

h — = dh~ ^v s ' ni 'W J de s S{e,e s ,a,fi) 

2 f 

J[ dej c J c (e Jc ,a,n)5(e-e Jl -ej 2 -e s ), (27) 



de de dh^ 

■1 

x 

c 
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with (as above) H = 1 + 0(a s ). Referring to the notation of Eqs. (21) and (22), the functions S 
and J c are: 

S{e,e 8 ,a,n) = ]T S(N s ,^)5(e- f(N s ))5(e s - f(N s ,a)) (28) 
J e (e Je ,a,n) = -(2nf ^J c (N Jc ^)S(e Jc - f c (N Jc: a))5 -u;(N Jc )) 5 2 (n x ± n(N Jc )), 

S N Jc \ A / 

(29) 

with the plus sign in the angular delta function for jet 2, and the minus for jet 1. The weight 
functions for the jets are given by Eq. (22) and induce dependence on the parameter a. We have 
introduced the factorization scale //, which we set equal to the renormalization scale. 

We note that we must construct the soft functions S(N s ,fi) to cancel the contributions of 
final-state particles from each of the J c (N Jc) //) to the weight e, as well as the contributions of 
the jet functions to e from soft radiation outside their respective regions Cl c . Similarly, the jet 
amplitudes must be constructed to include collinear enhancements only in their respective jet 
directions. Explicit constructions that satisfy these requirements will be specified in the following 
subsections. 

To disentangle the convolution in (27), we take Laplace moments with respect to e: 

da(e,u, s,a) r 00 , _„ = da(e,e,a) 

1 dee 



de dh\ Jo de de dfi\ 

da 



dhi 



H(s,h-y,n) S(e,v,a,ii) ][ J c (v, a, //). (30) 



c=l 



Here and below unbarred quantities are the transforms in e, and barred quantities denote un- 

transformed functions. _ 

f'°° _ — 

S(e,v, a,//) = / de s e~ ves S(e,e s ,a,n), (31) 
Jo 

and similarly for the jet functions. 

In the following subsections, we give explicit constructions for the functions participating in 
the factorization formula (27), which satisfy the requirement of infrared safety, and avoid double 
counting. An illustration of the cross section factorized into these functions is shown in Fig. 4. 
As discussed above, non-global logarithms will emerge when ev becomes small enough. 



3.3 The short-distance function 

The power counting described in [24] shows that in Feynman gauge the subdiagrams of Fig. 4 
that contribute to H in Eq. (27) at leading power in e and e are connected to each of the two jet 
subdiagrams by a single on-shell quark line, along with a possible set of on-shell, collinear gluon 
lines that carry scalar polarizations. The hard subdiagram is not connected directly to the soft 
subdiagram in any leading region. 
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Figure 4: Factorized cross section (27) after the application of Ward identities. The vertical line 
denotes the final state cut. 



The couplings of the scalar-polarized gluons that connect the jets with short-distance sub- 
diagrams may be simplified with the help of Ward identities (see, e.g. [26]). At each order of 
perturbation theory, the coupling of scalar-polarized gluons from either jet to the short-distance 
function is equivalent to their coupling to a path-ordered exponential of the gauge field, oriented 
in any direction that is not collinear to the jet. Corrections are infrared safe, and can be ab- 
sorbed into the short- distance function. Let h(pj c , hi, A) represent the set of all short-distance 
contributions to diagrams that couple any number of scalar-polarized gluons to the jets, in the 
amplitude for the production of any final state. The argument A stands for the fields that create 
the scalar-polarized gluons linking the short-distance function to the jets. On a diagram-by- 
diagram basis, h depends on the momentum of each of the scalar-polarized gluons. After the 
sum over all diagrams, however, we can make the replacement: 



where /12 is a short- distance function that depends only on the total momenta pj 1 and pj 2 . It 
also depends on vectors £ c that characterize the path-ordered exponentials $(0, —00; 0): 



where the superscript (f) indicates that the vector potential takes values in representation f, 
in our case the representation of a quark or antiquark. These operators will be associated with 
gauge-invariant definitions of the jet functions below. To avoid spurious collinear singularities, we 
choose the vectors £ c , c = 1, 2, off the light cone. In the full cross section (30) the ^-dependence 
cancels, of course. 




(32) 




(0, -00; 0) = Pe 



i9 fl^dX UA (l H^c) 



(33) 
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The dimensionless short- distance function H = \h 2 \ 2 in Eq. (27) depends on y/s and pj c ■ £ c , 
but not on any variable that vanishes with e and e: 

H(pj c ,£c,ni,n) = H (yH 7 PJc ^ c ,ni,a a (/x) ) , (34) 

where 

Ic^c/)/^. (35) 

Here we have observed that each diagram is independent of the overall scale of the eikonal vector 
3.4 The jet functions 

The jet functions and the soft functions in Eq. (27) can be defined in terms of specific matrix 
elements, which absorb the relevant contributions to leading regions in the cross section, and 
which are infrared safe. Their perturbative expansions specify the functions S and J c of Eq. 
(29). We begin with our definition of the jet functions. 

The jet functions, which absorb enhancements collinear to the two outgoing particles produced 
in the primary hard scattering, can be defined in terms of matrix elements in a manner reminiscent 
of parton distribution or decay functions [27]. To be specific, we consider the quark jet function: 

s Mc Njc 

x 5(e Jc - f c (N Jc , a)) 5 ^ - u(N Je ) j 5 2 (n c - n(N Jc )) , (36) 

where Mc is the number of colors, and where h c denotes the direction of the momentum of jet 
c, Eq. (29), with n 2 = —n\. q is the quark field, <3>^(0, — oo; 0) a path-ordered exponential in 
the notation of (33), and the trace is taken over color and Dirac indices. We have chosen the 
normalization so that the jet functions J' M in (36) are dimensionless and begin at lowest order 
with 

jy o \ej e ,a,n) = %6(ej e ), (37) 
with /3g the lightlike velocities corresponding to the jet momenta in Eq. (10): 

0f = Vn # = V- (38) 

The scalar jet functions of Eq. (29) are now obtained by projecting out the component of J'J 1 in 
the jet direction: 

J c (e Jc , a, n) = j3 c - J' c {£j c , a, n) = 5{e Jc ) + G(a s ) , (39) 

where (5\ — fa, fii — Pi are the lightlike vectors in the directions opposite to (5\ and fa, respec- 
tively. By construction, the J c are linear in fa. 



Y{0 $g }t (0,-oo;0M0) Nj)(Nj, 



g(0)^(0,-oo;0) 
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To resum the jet functions in the variables Ej c , it is convenient to reexpress the weight 
functions (22) in combinations of light-cone momentum components that are invariant under 
boosts in the X3 direction, 

h(N j2 ,a) = -±- 2 J2 K± ( 2 P~jXY~ a - (40) 

Here we have used the relation \/s/2 = u Jc , valid for both jets in the cm. At the same time, we 
make the identification, 

-J - cv(N Jc )j 5\n c - h{N Jc )) = \ 5 3 (p Jc - p(N Jc )) , (41) 

which again holds in the cm. frame. The spatial components of each pj c are thus fixed. Given that 
we are at small e,j c) the jet functions may be thought of as functions of the light- like jet momenta 
Pj c of Eq. (10) and of e j c . Because the vector jet function is constructed to be dimensionless, 
in Eq. (36) is proportional to f3 c rather than pj c . Otherwise, it is free of explicit j3 c - dependence. 

The jet functions can now be written in terms of boost-invariant arguments, homogeneous of 
degree zero in £ c : 



Jc [£jci a i A 4 ) — Pen 



% J c (1) 



ntA^c 7(2) 
led 2 c 




(42) 



where and are independent functions, and where we have suppressed possible dependence 
on £ c ,j_. For jet c, the weight e,j c is fixed by S(e j c — f c (Nj c , a)), where on the right-hand side of 
the expression for the weight (40), the sum over each particle's momentum involves the overall 
factor (2pj c / v /i) 1 ~ a . After integration over final states at fixed ej c) the jet can thus depend on 
the vector pj c . At the same time, it is easy to see from the definition of the weight that pj c can 
only appear in the combination (l/£ r j c - v /s) 1 /( 1_a ) (2p% Jy/s). This vector can combine with £ c to 
form an invariant, and all ^-dependence comes about in this way. 
Expression (42) can be further simplified by noting that 

2&-£ c /W c = C 2 + C 2 ± • (43) 
Choosing £ Cj j_ = 0, we find a single combination, 

J c (e Jc , a, /i) = J c ( *±Jz, e Jc ^ (Ce) 1 ^ , a, , (44) 
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where, in the notation of Eq. (42), J c = + j( 2 \ and we have defined 

In these terms, the Laplace moments of the jet function inherit dependence on the moment 
variable v through 



r°° _ _ 

Jc (v,a,n) = / de Jc e~ U£j - J c (e Jc , a, //) 
Jo 



= J C { P -^A (Cc) 1_a , a, a.(|i)V (46) 

where the unbarred and barred quantities denote transformed and untransformed functions, 
respectively. We have constructed the jet functions to be independent of e, since the radiation 
into Q is at wide angles from the jet axes and can therefore be completely factored from the 
collinear radiation. This radiation at wide angles is contained in the soft function, which will be 
defined below in a manner that avoids double counting in the cross section. 



3.5 The soft function 



Given the definitions for the jet functions in the previous subsection, and the factorization (27), 
we may in principle calculate the soft function S order by order in perturbation theory. We can 
derive a more explicit definition of the soft function, however, by relating it to an eikonal analog 
ofEq. (27). 

As reviewed in Refs. [15, 26], soft radiation at wide angles from the jets decouples from the 
collinear lines within the jet. As a result, to compute amplitudes for wide-angle radiation, the jets 
may be replaced by nonabelian phases, or Wilson lines. We therefore construct a dimensionless 
quantity, a^ eik \ in which gluons are radiated by path-ordered exponentials $, which mimic the 
color flow of outgoing quarks, 



(47) 



with p c a light-like velocity in either of the jet directions. For the two-jet cross section at measured 
e and e cik , we define 



(e, £ eik , a, = _ £ (o $g jt (oo, 0; 0)Q%*(oo, 0; 0) 



x N 



(oo, 0; 0)$g>(oo, 0; 0)| 0) 5 (e - /(7V eik )) 5 (e cik - f(N ci 



6(e) S(e eik ) + 0(a s ) . 



(48) 



The sum is over all final states A^ eik in the eikonal cross section. The renormalization scale in 
this cross section, which will also serve as a factorization scale, is denoted ji. Here the event 
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shape function e cik is denned by f(N eik ,a) as in Eqs. (3) and (4), distinguishing between the 
hemispheres around the jets. As usual, Mc is the number of colors, and a trace over color is 
understood. 

The eikonal cross section (48) models the soft radiation away from the jets, including the 
radiation into Q, accurately. It also contains enhancements for configurations collinear to the 
jets, which, however, are already taken into account by the partonic jet functions in (27). Indeed, 
(48) does not reproduce the partonic cross section accurately for collinear radiation. It is also 
easy to verify at lowest order that even at fixed e cik the eikonal cross section (48) is ultraviolet 
divergent in dimensional regularization, unless we also impose a cutoff on the energy of real gluon 
emission collinear to f3\ or (3%. 

The construction of the soft function S from is nevertheless possible because the eikonal 
cross section (48) factorizes in the same manner as the cross section itself, into eikonal jet functions 
and a soft function. The essential point [14] is that the soft function in the factorized eikonal 
cross section is the same as in the original cross section (27). The eikonal jets organize all 
collinear enhancements in (48), including the spurious ultraviolet divergences. These eikonal jet 
functions are defined analogously to their partonic counterparts, Eq. (36), but now with ordered 
exponentials replacing the quark fields, 

jt ik) (ec,a,ri = -i- £ (0 Isg^O, -oo; 0)^(00,0; 0)1 

( N W ^^(oo, 0; 0)<t>£ } (0, -00; 0)\0)6 (e c - f c (N^\a)) 
= 5(e c ) + 0(a s ), (49) 

where f c is a quark or antiquark, and where the trace over color is understood. The weight 
functions are given as above, by Eq. (22), with the sum over particles in all directions. 
In terms of the eikonal jets, the eikonal cross section (48) factorizes as 

2 

a- (oik) (e, e cik , a, fi) = J de s S (e, e 8 , a,/j,)f[ J de c J c (oik) (e c , a, fi) 5 (e cik - e s - e x - e 2 ) , (50) 

c=l 

where we pick the factorization scale equal to the renormalization scale \x. As for the full cross 
section, the convolution in (50) is simplified by a Laplace transformation (46) with respect to 
e cik , which allows us to solve for the soft function as 

S (e, 1/, a, pt) = f ^ (£ ' a ' ^ = 5(e) + 0(a s ) . (51) 
FT. Jc (oik) (1/, a, 

c=l 

In this ratio, collinear logarithms in v and the unphysical ultraviolet divergences and their associ- 
ated cutoff dependence cancel between the eikonal cross section and the eikonal jets, leaving a soft 
function that is entirely free of collinear enhancements. The soft function retains //-dependence 
through soft emission, which is also restricted by the weight function e. In addition, because 
soft radiation within the eikonal jets can be factored from its collinear radiation, just as in the 
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partonic jets, all logarithms in v associated with wide-angle radiation are identical between the 
partonic and eikonal jets, and factor from logarithmic corrections associated with collinear radi- 
ation in both cases. As a result, the inverse eikonal jet functions cancel contributions from the 
wide-angle soft radiation of the partonic jets in the transformed cross section (30). 

Given the definition of the energy flow weight function /, Eq. (5), the soft function is not 
boost invariant. In addition, because it is free of collinear logs, it can have at most a single 
logarithm per loop. Its dependence on e is therefore only through ratios of the dimensional 
quantities s^Js with the renormalization (factorization) scale. 

As in the case of the partonic jets, Eq. (46), we need to identify the variable through which 
v appears in the soft function. We note that dependence on the velocity vectors f3 c and the 
factorization vectors £ c must be scale invariant in each, since they arise only from eikonal lines 
and vertices. The eikonal jet functions cannot depend explicitly on the scale-less, lightlike eikonal 
velocities f3 C) and is independent of the £ c . Dependence on the factorization vectors £ c enters 
only through the weight functions, (40) for the eikonal jets, in a manner analogous to the case of 
the partonic jets. This results in a dependence on (Cc) 1_a , as above, with ( c defined in Eq. (45). 
In summary, we may characterize the arguments of the soft function in transform space as 

S (e, v, a,/i) = S ( — , ev, — (Q l ~ a , a, a s (/i) j . (52) 



4 Resummation 

We may summarize the results of the previous section by rewriting the transform of the factorized 
cross section (30) in terms of the hard, jet and soft functions identified above, which depend on 
the kinematic variables and the moment v according to Eqs. (34), (46) and (52) respectively, 

da(e,v,s,a) da ( ' y/s p Jc ■ | c . \ * (pj c -L \fs x _ a \ 

dedni drii \ [i /j, J " y /x [iv J 

xs(^,ev,£(ti 1 -\a,a.(jjL)) . (53) 

The natural scale for the strong coupling in the short-distance function H is \fsj1. Setting 
\i = y/s/2, however, introduces large logarithms of e in the soft function and large logarithms of 
v in both the soft and jet functions. The purpose of this section is to control these logarithms 
by the identification and solution of renormalization group and evolution equations. 

The information necessary to perform the resummations is already present in the factorization 
(53). The cross section itself is independent of the factorization scale 

d da(eus,a) = 
dfi dedrii 

and of the choice of the eikonal directions, £ c , used in the factorization, 

d da (e, v, s, a) 



d\n(p Jc -i c ) dedn x 



. (55) 
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The arguments of this section closely follow the analysis of Ref. [30]. We will see that the 
dependence of jet and soft functions on the parameter a that characterizes the event shapes 
(3) is reflected in the resummed correlations, so that the relationship between correlations with 
different values of a is both calculable and nontrivial. 



4.1 Energy flow 

As a first step, we use the renormalization group equation (54) to organize dependence on the 
energy flow variable e. Applying Eq. (54) to the factorized correlation (53), we derive the following 
consistency conditions, which are themselves renormalization group equations: 

H^j- In S ( ^-,ev, — (( c y~ a ,a,a s {fi) J = -7, (a a (fi)) , (56) 
d(JL \ \x \xv ) 

% ^ Jc (^^'^ (Cc)l " a ' a ' as(/i) ) = -^e(«.M). (57) 

ln H \ —, PJc — ,wi,a a (^)) = 7 s (a s (//)) + J2"fJc( a s(n)) ■ (58) 

The anomalous dimensions 7^, d = s, J c can depend only on variables held in common between 
at least two of the functions. Because each function is infrared safe, while ultraviolet divergences 
are present only in virtual diagrams, the anomalous dimensions cannot depend on the parameters 
v, e or a. This leaves as arguments of the 7^ only the coupling a s (/i), which we exhibit, and ( c , 
which we suppress for now. 

Solving Eqs. (56) and (57) we find 



S I -*—,ev, — (Q L a ,a,a s (n)\ = Si ,ev, (C c ) a ,a,a s (/j, ) e "° 



(59) 



Jc[ , — (Cc) ,a,a a (n)\ = Jc\— , — (Cc) ,a,a a (fio) e m 



J f 7Jc ( Qs (A)) 



(60) 



for the soft and jet functions. As suggested above, we will eventually pick ji ~ y/s to avoid large 
logs in H. Using these expressions in Eq. (53) we can avoid logarithms of e or v in the soft 
function, by evolving from /x = Sy/s to the factorization scale ji ~ y/s. No choice of /2 , however, 
controls all logarithms of v in the jet functions. Leaving p, free, we find for the cross section 
(53) the intermediate result 

da (e, v, s, a) da ( y/s p Jc ■ | c A \ 
dedrii dn\ \ [i fx I 
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X 



S (l,eu, (Cc) 1 a ,a,a s (ey/s)^j exp 




dX 

T 



7 s (a s (A)) > 



(61) 



l-a 



,a,a s (p, ) exp 




<2A 
T 



7j c (a«(A)) > ■ 



We have avoided introducing logarithms of e into the jet functions, which originally only depend 
on u, by evolving the soft and the jet functions independently. The choice of /io = £\/s or \fsjv 
for the soft function is to some extent a matter of convenience, since the two choices differ by 
logarithms of ev. In general, if we choose /j,q = y/s/v, logarithms of ev will appear multiplied by 
coefficients that reflect the size of region f2. An example is Eq. (15) above. When f2 has a small 
angular size, fiQ = \fsjv is generally the more natural choice, since then logarithms in ev will 
enter with small weights. In contrast, when Q grows to cover most angular directions, as in the 
study of rapidity gaps [32], it is more natural to choose /i = £\/s. 

4.2 Event shape transform 

The remaining unorganized "large" logarithms in (61), are in the jet functions, which we will 
resum by using the consistency equation (55). The requirement that the cross section be in- 
dependent of pj c ■ £ c implies that the jet, soft and hard functions obey equations analogous to 
(56)-(58), again in terms of the variables that they hold in common [30]. The same results may 
be derived following the method of Collins and Soper [21], by defining the jets in an axial gauge, 
and then studying their variations under boosts. 

For our purposes, only the equation satisfied by the jet functions [21, 30] is necessary, 



The functions K c and G c compensate the ^-dependence of the soft and hard functions, respec- 
tively, which determines the kinematic variables upon which they may depend. In particular, 
notice the combination of v- and ^-dependence required by the arguments of the jet function, 



Since the definition of our jet functions (36) is gauge invariant, we can derive the kernels 
K c and G c by an explicit computation of dJ c /d\n [p,j c ■ £ c ) in any gauge. The multiplicative 
renormalizability of the jet function, Eq. (57), with an anomalous dimension that is independent 
of pj c -£ c ensures that the right-hand side of Eq. (62) is a renormalization-group invariant. Thus, 
K c + G c are renormalized additively, and satisfy [21] 




(62) 



Eq. (46). 
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an 



(63) 



Since G c and hence 7^ c , may be computed from virtual diagrams, they do not depend on a, and 
7x c is the universal Sudakov anomalous dimension [21, 33]. 

With the help of these evolution equations, the terms K c and G c in Eq. (62) can be reexpressed 
as [34] 



K c [ (Cc) 1 a , a, a s (p)j + G c f PJc ^ > ««(a0 



K c [—,a,a s I ci 



^(Q'-^+Gc (la, (c 2 p Jc -cf c ) 



C2 PJ C 'Cc 



^-7* e («.(*')) 



C2PJ c -£c 



= -B' c (c 1 ,c 2 ,a,a a (c 2 pj c -Q)-2 J —A' c (c 1 ,a,a s (X')) , 



(64) 



ciV^CCc) 1 " /" 



where in the second equality we have shifted the argument of the running coupling in K c , and 
have introduced the notation 



1 



B' c (c 1 ,c 2 ,a,a s {fi)) = -K c [—,a,a s {fi))-G c [—,a s {fi)), 



Cl 



c 2 



2^(ci,a,a s (//)) = 7^ c (a s (//)) + P(g(fJ,))-^-^K c (— ,a,a s (//)J . 



(65) 



The primes on the functions A' c and 1^ are to distinguish these anomalous dimensions from their 
somewhat more familiar versions given below. 
The solution to Eq. (62) with /i = p, is 



Jc 



'VJ C -ic Vs u _ a 



(Cc) a ,a,a s (n )\ = J c 



V* V* /A \l-a 



x exp 



v^/(2Co) 



T 



v 2 Co Mo 

C2 A 



(Co) a ,a,a s (/i ) 
^' A , 



/Cla 
~Y A 'c {ci,a,a s (A')) 



1 ,(2A)!-» 



(66) 



where we evolve from y/s/ (2 Co) to pj c ■ £ c = yfs/ (2 Cc) (see Eq. (45)) with 

, x l/(2-o) 

Co ' ' ' 



(67) 



After combining Eqs. (60) and (66), the choice /x = y/s/(2( ) = ^(Co) 1 a allows us to control 
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all large logarithms in the jet functions simultaneously: 



. 1 



Jr 



(Cc) a ,a,a s (/j,)\ = J c [l,l,a,a a 



2 Co, 



exp < 



dX 
T 



7j c MX)) 



I v^/(2Co) 



x exp < — 



dX 
T 



C2 A 



/■ dA' 

5c(ci,cj,o,a ( (cjA)) + 2 y — Aj.(ci,a,a a (A')) 
As observed above, we treat a as a fixed parameter, with |a| small compared to In (1/e) and lnz/. 



(68) 



4.3 The resummed correlation 

Using Eq. (68) in (61), and setting jx = y/s/2, we find a fully resummed form for the correlation, 

da (e, v, s, a) da ^ ( 2p Jc -j c ( \T S X\ 

dedfii dhi \ y/s ' ' s \ 2 J J 



x S(l,ev, (Q 1 - a ,a,a 8 (£y/s)) exp I - J —-? s (a s (X)) 



x 



c=l V 



2 Co, 



exp < 



v^/2 



I v^/(2Co) 



/ yU(« s (A)) 



PJ c <c 



x exp < — 



dX 

T 



v^/(2Co) 



C /- A dA' 

B' c (c 1 ,c 2 ,a,a s (c 2 X)) + 2 y — A' c (ci, a, a s (A')) 



1 „(2A)l-« 



(69) 



Alternatively, we can combine all jet-related exponents in Eq. (69) in the correlation. As we 
will verify below in Section 5.2, the cross section is independent of the choice of £ c . As a result, 
we can choose 

PJ C ■ i = ^ ■ (70) 
This choice allows us to combine 7j c and B' c in Eq. (69), 

da (g, z/, s, a) d<7o / / 

de dni dni \ ' \ 2 / / 

1 After this paper was submitted for publication, a related analysis of event shape and energy flow correlations 
was given by Dokshitzer and Marchesini [31], who identify the same factorization of soft radiation described here 
and in [18], and who study the leading logarithms of ev for e <C l/v, using the methods of [20]. 
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x S (l, eu, 1, a, a s (e\^s)^j exp 

v^/2 



y/5/2 



dX 

T 



7«(«-W) 



x exp < 



T 



v^/(2Co) 



C2 / dA' 

7J c (a s (A)) + B' c (c 1 ,c 2 ,a,a s (c 2 \)) + 2 y — (c u a, a s (A')) 



1 !/(2 A)l-" 



(71) 



with ( Q given by Eq. (67). 

In Eqs. (69) and (71), the energy flow e appears at the level of one logarithm per loop, in S, 
in the first exponent. Leading logarithms of e are therefore resummed by knowledge of 7^, the 
one-loop soft anomalous dimension, where we employ the standard notation, 



7* (« 



.) = Et 



(n) 



n=0 



(72) 



for any expansion in a s . At the same time, v appears in up to two logarithms per loop, char- 
acteristic of conventional Sudakov resummation. To control //-dependence at the same level as 
^-dependence, it is natural to work to next-to-leading logarithm in u, by which we mean the 
level a k s ln fc v in the exponent. As usual, this requires one loop in B' c and 7j c , and two loops in 
the Sudakov anomalous dimension A' c , Eq. (65). These functions are straightforward to calculate 
from their definitions given in the previous sections. Only the soft function S in Eqs. (69) and 
(71) contains information on the geometry of Q. The exponents are partially process-dependent, 
but geometry- independent. In Section 5, we will derive explicit expressions for these quantities, 
suitable for resummation to leading logarithm in e and next-to-leading logarithm in v. 



4.4 The inclusive event shape 

It is also of interest to consider the cross section for e + e~-annihilation into two jets without fixing 
the energy of radiation into Q, but with the final state radiation into all of phase space weighted 
according to Eq. (4), schematically 



e 



+ 



+ e -> Jiipj^ + J 2 (pj 2 , /n 2 ) , (73) 



where Cli and Cl 2 cover the entire sphere. This cross section can be factorized and resummed in 
a completely analogous manner. The final state is a convolution in the contributions of the jet 
and soft functions to e as in Eq. (27), but with no separate restriction on energy flow into fl 
All particles contribute to the event shape. We obtain an expression very analogous to Eq. (69) 
for this inclusive event shape in transform space, which can be written in terms of the same jet 
functions as before, and a new function S mcl for soft radiation as: 

da incl (u, s, a) da (2p Jc ■ | c A \ 
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x S""-' ( (c,,.) 1 " (^j j exp<|- y y 7s (« s (A)) 

v/i/j/ 



dA 



x 



c=1 V V /y I ^i/(2 Co) J 



x exp < 



PJ c <c 

- / 

V5/(2Co) 



dA 

T 



c f dX' 

B' c (ci,C2,a, a s (c 2 A)) + 2 y — A' c (ci, a, a s (A')) 



1 >/(2 A) 1 "" 



(74) 



Here the soft function S mcl = 1 + 0(a s ). The double-logarithmic dependence of the shape 
transform is identical to our resummed correlation, Eq. (69). We will show below, in Sec. 5.3, 
that Eq. (74) coincides at NLL with the known result for the thrust [6] when we choose a = 0. 

5 Results at NLL 

5.1 Lowest order functions and anomalous dimensions 

In this section, we describe the low-order calculations and results that provide explicit expressions 
for the resummed shape/flow correlations and inclusive event shape distributions at next-to- 
leading logarithm in v and leading logarithm in e (we refer to this level collectively as NLL 
below). We go on to verify that for the case a = we rederive the known result for the resummed 
thrust at NLL, and we exhibit the expressions for the correlation that we will evaluate in Sec. 6. 



5.1.1 The soft function 

The one-loop soft anomalous dimension is readily calculated in Feynman gauge from the combi- 
nation of virtual diagrams in <r( cik ), Eq. (48), and J (eik) , Eq. (49), in Eq. (51). The calculation 
and the result are equivalent to those of Ref . [14] , where the soft function was formulated in axial 
gauge, 



£ In (ft • i) - In 



- 1 



c=l 



(75) 



The first, ^-dependent logarithmic term is associated with the eikonal jets, while the second 
is a finite remainder from the combination of <r( clk ) and j( eik ) in (51). Whenever £ C) j_ = 0, the 
logarithmic terms cancel identically, leaving only the final term, which comes from the £ c eikonal 
self-energy diagrams in the eikonal jet functions. 

The soft function is normalized to S^°\e) = 5(e) as can be seen from (51). For non-zero e, 
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da/de is given at lowest order by 

S w (e ± 0, n) = C F - ( dPS 2 P\'^ 2 , , (76) 

where PS2 denotes the two-dimensional angular phase space to be integrated over region Q, and 
k = k/u k . We emphasize again that the soft function contains the only geometry-dependence of 
the cross section. Also, for e 7^ is independent of v and a. 

As an example, consider a cone with opening angle 25, centered at angle a from jet 1. In this 
case, the lowest-order soft function is given by 

5« (e ? 0, a, S) = C F l - In f . (77) 

e \ cos z a — cos^ J 

Similarly, we may choose Q as a ring extending angle Si to the right and 5 2 to the left of the 
plane perpendicular to the jet directions in the center-of-mass. In this case, we obtain 

r.m / /n r r\ ^ 1, / (1 + sin 5i) (1 + sin 5 2 ) \ ^ 2 

e ^ 0, S u S 2 ) = C F - In )- — {)- — { = C F - A?? , 78 

£ \(1 — smOi) (1 — sin^j y e 

with A77 the rapidity spanned by the ring. For a ring centered around the center-of-mass (Si = 
S2 = S) the angular integral reduces to the form that we encountered in the example of Sec. 2.3, 
and that we will use in our numerical examples of Sec. 6, with Ar] given by Eq. (14). 

5.1.2 The jet functions 

Recall from Eq. (39) that the lowest-order jet function is given by jjP^ = 1. 
The anomalous dimensions of the jet functions are found to be 

-y^ = ~c F , (79) 

the same for each of the jets. The jet anomalous dimensions are process-independent, but of 
course flavor-dependent. The same anomalous dimensions for final-state quark jets appear in 
three- and higher-jet cross sections. 

5.1.3 The A-G-decomposition 

The anomalous dimension for the A-G-decomposition is, as noted above, the Sudakov anomalous 
dimension, 

7g = 2 C F , (80) 
7g = KC F , (81) 
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also independent of the jet-direction. The well-known coefficient K (not to be confused with the 
functions K c ) is given by [35] 

with the normalization T F — 1/2 and Nj the number of quark flavors. 

K c and G c , the functions that describe the evolution of the jet functions in Eq. (62), are given 
at one loop by 



(83) 



Evolving them to the values of pL with which they appear in the functions A' c and B' c , Eq. (65), 
they become 

= ~Cf In (e^-^c?) , (85) 
= -^Infe-lV (86) 



c 2 J V r : 



Recall that G c is computed from virtual diagrams only, and thus does not depend on the weight 
function. It therefore agrees with the result found in [21]. The soft-gluon contribution, K c , which 
involves real gluon diagrams, does depend on the cross section being resummed. 
With the definitions (65) of A' c and B' c we obtain 



A'V = C F , (87) 
K+^\n (e 2 ^- 1+a cj] 
2d 



A>W(c 1: a) = \c F 



(88) 

B>V( Cl ,c 2 ,a) = 2C F \n(e^- 1+a ^y (89) 



Here (3q is the one- loop coefficient of the QCD beta-function, A) = 5 (ILA/c - 4T F N f ) ((3(g) = 
-9tPo + 0(g 3 )). 

5.1.4 The hard scattering, and the Born cross section 

At NLL only the lowest-order hard scattering function contributes, which is normalized to 

ff(°>(a,(Vi/2)) = l. (90) 
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At this order the hard function is independent of the eikonal vectors £ c , although it acquires 
^-dependence at higher order through the factorization described in Sec. 3.3. For completeness, 
we also give the electromagnetic Born cross section at fixed polar and azimuthal angle: 

t=M?«)£( 1+ «"•')• (91) 

where 6 is the cm. polar angle of hi, eQf is the charge of quark flavor f, and a crn = e 2 /(47r) is 
the fine structure constant. 



5.2 Checking the ^-dependence 



It is instructive to verify how dependence on the eikonal vectors £ c cancels in the exponents of 
the resummed cross section (69) at the accuracy at which we work, single logarithms of e, and 
single and double logarithms of v. In these exponents, ^-dependence enters only through the 
combinations (f3 c ■ £ c ) and (pj c ■ £ c ). 

Let us introduce the following notation for the exponents in Eq. (69), to which we will return 
below: 



^5/2 2 V^/2 

Ei = ~ J y7. («.(A)) - E / y7J e («.(A)), 



S\fs 

2 PJci 



v^/(2Co) 



E 

c=l 



dX 
T 



v^/(2Co) 



ci A 

5^ (ci,c 2 ,a,a s (c 2 A)) + 2 / 



Cl 



k(2 A) 1 " 



dA' 

— A' c (ci,a,a s (X')) 



(92) 



(93) 



At NLL, explicit £ c dependence is found only in 7 S , Eq. (75), for Ei, and in the upper limit of 
the A integral of E 2 . We then find that 



d 



dln/3 c .£ 



(Ei + E 2 ) = 2C F 



£V ^ 1 ,(2 PJc .C c )l-« 



Here the second term stems entirely from A'^\ Eq. (87); other contributions of E 2 are sub- 
leading. The ^-dependence in the exponents begins only at the level that we do not resum, 
at a s ln(l/ez/), which is compensated by corrections in S(eu,a s ). The remaining contributions 
are of NNLL order, that is, proportional to at k s (^fs) ln fc_1 (u (3 C -£ c ), as may be verified by ex- 
panding the running couplings. Thus, as required by the factorization procedure, the relevant 
^-dependence cancels between the resummed soft and jet functions, which give rise to the first 
and second integrals, respectively, in Eq. (94). 
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5.3 The inclusive event shape at NLL 

We can simplify the differential event shape, Eq. (74), by absorbing the soft anomalous dimension 
7 S into the remaining terms. We will find a form that can be compared directly to the classic 
NLL resummation for the thrust (a = 0). This is done by rewriting the integral over the soft 
anomalous dimension as 

\ r, T 7sQ!sA = r)I , M1/( , ,i X A + / ni -Yls(a s (X)) 

J^/s/u A J v /s/[2(i//2) 1 /(2-a)j \ J^/s/v A 

= ( 2 -°) frt, /«/« )1 ~Y 7s( a s(X)) 

J^/l/[2(u/2) 1 /( 2 - a )] A 

rVs/2 d\ r x d\' d 

In the first equality we split the A integral so that the limits of the first term match those of 
the B' c integral of Eq. (74). In the second equality we have changed variables in the second term 
according to 

/ s i-§ \ — 

so that the limits of the second integral also match. In the third equality of Eq. (95), we 
have reexpressed the running coupling at the old scale A in terms of the new scale. This is a 
generalization of the procedure of Ref . [36] , applied originally to the threshold- resummed Drell- 
Yan cross section [37]. 

Using Eq. (95), and identifying p Jc ■ £ c with y/s/2 (Eq. (70)) in the inclusive event shape 
distribution, Eq. (74), we can rewrite this distribution at NLL as 

da incl (z/, s, a 



x d\' 




/dX 
— A c (d,a,a s (A')) 



1 „(2A)l-« 



(97) 

where we have rearranged the contribution of 7 S as: 

1 d 

A c (c u a,a s (fi)) = A' c (c 1 ,a,a s (fi))--(l-a)P(g(ti))—'y a (a 8 (fi)), 

B c (c 1 ,c 2 ,a,a s (fj,)) = 7j c (a 8 (fj,)) + (l - |) 7 S (a 8 (fj,)) + B' c (a, c 2 , a, a s (//)) . (98) 
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Next, we replace the lower limit of the A'-integral by an explicit ^-function. Then we exchange 
orders of integration, and change variables in the term containing A from the dimensionful 
variable A to the dimensionless combination 



2AA' 



We find 

da incl (u, s, a) 
dhi 



da -X 



x ]~J exp 

c=l 



U 



(99) 



T 



B c (d,c 2 ,a,a s (A)) 



2(z//2)V(2- ) 



A' A' 2 /° 7/ v 1 u 



la„,\—a 



U 



/2 



lj A c ( Cl ,a,a s (A'))j . 

(100) 



Here, the ^-function vanishes for small A', and the remaining effects of replacing the lower bound- 
ary of the A' integral by are next-to- next-to-leading logarithmic. 

A further change of variables allows us to write the NLL resummed event shapes in a form 
familiar from the NLL resummed thrust. In the first line of Eq. (100), we replace A 2 — > us/ '4. In 
the second line we relabel A' — > y/q 2 , and exchange orders of integration. Finally, choosing 



we find at NLL 



da mcl (u, s, a) 
dh\ 



da 2 
dm 



Cl 

c 2 



i r us 
du 



-7b 



(101) 



+ i5 c (« s (W4))(e-^ 2/(2 ^— -1 



and reproduce the well-known coefficients 



Cf, 

^ = \C F K, 



3 c 



independent of a. In Eq. (102), we have made use of the relation 

e -*/!/ _ i « -0 (x-ye~ lE ) , 



(102) 

(103) 
(104) 

(105) 
(106) 



31 



which is valid at NLL in the logarithmic integrals. With these choices, when a = we reproduce 
the NLL resummed thrust cross section [6]. 

The choices of the q in Eq. (101) cancel all purely soft NLL components (7 S and K c ). The 
remaining double logarithms stem from simultaneously soft and collinear radiation, and single 
logarithms arise from collinear configurations only. At NLL, the cross section is determined 
by the anomalous dimension A c , which is the coefficient of the singular 1/[1 — x] + term in the 
nonsinglet evolution kernel [38], and the quark anomalous dimension. All radiation in dijet events 
thus appears to be emitted coherently by the two jets [6]. This, however, is not necessarily true 
beyond next-to-leading logarithmic accuracy for dijets, and is certainly not the case for multijet 
events [14]. Similar considerations apply to the resummed correlation, Eq. (69). 



5.4 Closed expressions 



Given the explicit results above, the integrals in the exponents of the resummed correlation, Eq. 
(69), may be easily performed in closed form. We give the analytic results for the exponents of 
Eq. (69), as defined in Eqs. (92) and (93). As in Eq. (70), we identify pj c ■ £ c with \fsj2. 



4C F 

'a a (y/s/2)\ * I '»> 



a 



6Cp_ 

00 



a 



(107) 



1 4C F 



a. 



( £2V^\ I 

V 2Co J J 



a, 



a, 



( ci Vs \ I 
V 2Co J / 



(108) 



with 



Ki(a) = In 



4 \ 4tt 



47T 



. 2Co y 



-i -1 



2K A 



In 



K 2 (a) = (1 - a - 2>y E ) + — 

Po 



a.. 



n -l 



Po 2(31 



2K (3i 



'A' 
Aire 



C 2 \fs\ ( C 2 y/s'' 



(109) 



In 



A) 



'ci 



ci v« 



2Co 



(110) 



We have used the two-loop running coupling, when appropriate, to derive Eqs. (107) - (110). 
The results are expressed in terms of the one-loop running coupling 

2tt 1 



and the first two coefficients in the expansion of the QCD beta-function, (3o and 

Pi = Y<%-(y Ca + 4Cf ) TpNf - 



(111) 



(112) 
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Combining the expressions for the exponents, Eqs. (107) and (108), for the Born cross section, 
Eq. (91), and for the soft function, Eq. (76), in Eq. (69), the complete differential cross section, 
at LL in e and at NLL in u, is given by 

da(e,u,s,a) f \ na 2 cin f 2 \ a 8 (ey/s) 1 r _ c 1 A • (3 2 



6C F 





These are the expressions that we will evaluate in the next section. We note that this is not the 
only possible closed form for the resummed correlation at this level of accuracy. When a full 
next-to-leading order calculation for this set of event shapes is given, the matching procedure of 
[6] may be more convenient. 



6 Numerical Results 

Here we show some representative examples of numerical results for the correlation, Eq. (113). 
We pick the constants q as in Eq. (101), unless stated otherwise. The effect of different choices 
is nonleading, and is numerically small, as we will see below. In the following we choose the 
region Q to be a ring between the jets, centered in their center-of-mass, with a width of Arj = 2, 
or equivalently, opening angle 5 ~ 50 degrees (see Eq. (14)). The analogous cross section for a 
cone centered at 90 degrees from the jets (Eq. (77)) has a similar behavior. In the following, the 
center-of-mass energy Q = y/s is chosen to be 100 GeV. 

Fig. 5 shows the dependence of the differential cross section (69), multiplied by e and nor- 
malized by the Born cross section, £f ^^" 1 - ) , on the measured energy e and on the parameter a, 

at fixed v. In Fig. 5 a), we plot for v = 10, in Fig. 5 b) for v = 50. As v increases, the 

radiation into the complementary region is more restricted, as illustrated by the comparison 
of Figs. 5 a) and b). Similarly, as a approaches 1, the cross section falls, because the jets are 
restricted to be very narrow. On the other hand, as a assumes more and more negative values at 
fixed e, the correlations (69) approach a constant value. For a large and negative, however, non- 
global dependence on Ine and \a\ will emerge from higher order corrections in the soft function, 
which we do not include in Eq. (113). 

In Fig. 6 we investigate the sensitivity of the resummed correlation, Eq. (113), to our choice 
of the constants q. The effect of these constants is of next-to- next-to-leading logarithmic order 
in the event shape. We plot the differential cross section e ed ^^\ at Q = 100 GeV, for fixed 
e = 0.05 and v = 20, as a function of a. The effects of changes in the q are of the order of a few 
percent for moderate values of a. 
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o.i- 1 - 5 



Figure 5: Differential cross section £ J^y^" 1 , normalized by the Born cross section, at Q = 100 
GeV, as a function of £ and a at fixed v: a) z/ = 10, b) z/ = 50. Q is a ring (slice) centered around 
the jets, with a width of At] = 2. 



Finally, we illustrate the sensitivity of these results to the flavor of the primary partons. 
For this purpose we study the corresponding ratio of the shape/flow correlation to the cross 
section for gluon jets produced by a hypothetical color singlet source. Fig. 7 displays the ratio of 
the differential cross section da q (e,a)/(dedni), Eq. (113), normalized by the lowest-order cross 
section, to the analogous quantity with gluons as primary partons in the outgoing jets, again at 
Q = 100 GeV. This ratio is multiplied by Ca/Cf in the figure to compensate for the difference in 
the normalizations of the lowest-order soft functions. Gluon jets have wider angular extent, and 
hence are suppressed relative to quark jets with increasing v or a, as can be seen by comparing 
Figs. 7 a) and b). Fig. 7 a) shows the ratio at v — 10, and Fig. 7 b) at v — 50. These results 
suggest sensitivity to the more complex color and flavor flow characteristic of hadronic scattering 
[14,15]. 



7 Summary and Outlook 

We have introduced a general class of inclusive event shapes in e + e~ dijet events which reduce 
to the thrust and the jet broadening distributions as special cases. We have derived analytic 
expressions in transform space, and have shown the equivalence of our formalism at NLL with 
the well-known result for the thrust [6]. Separate studies of this class of event shapes in the 
untransformed space, at higher orders, and for nonperturbative effects [9] are certainly of interest. 
We reserve these studies for future work. 

We have introduced a set of correlations of interjet energy flow for the general class of event 
shapes, and have shown that for these quantities it is possible to control the influence of sec- 
ondary radiation and nonglobal logarithms. These correlations are sensitive mainly to radiation 
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Figure 6: Differential cross section eda J^^ ^ ' norma hzed by the Born cross section, at Q = 100 
GeV, as a function of a at fixed v = 20 and e = 0.05. Q is chosen as in Fig. 5. Solid line: 
ci = e~ 7E , C2 = 2, as in Eq. (101), dashed line: c\ = C2 = 1, dotted line: c\ = C2 = 2. 



emitted directly from the primary hard scattering, through transforms in the weight functions 
that suppress secondary, or non-global, radiation. We have presented analytic and numerical 
studies of these shape/flow correlations at leading logarithmic order in the flow variable and at 
next-to-leading-logarithmic order in the event shape. The application of our formalism to mul- 
tijet events and to scattering with initial state hadrons is certainly possible, and may shed light 
on the relationship between color and energy flow in hard scattering processes with non-trivial 
color exchange. 
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A Eikonal Example 

In this appendix, we give details of the calculation of the logarithmic behavior in the diagrams 
of Fig. 3. We choose the reference frame such that the momenta of the final state particles are 
given by: 

/3i = (1,0,0,1), 
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Figure 7: Ratios of differential cross sections for quark to gluon jets 

% ( ^Vdftf ) ( ^dJyS^ y 1 at Q = 100 GeV as a function of e and a at fixed i/: a) 
v — 10, b) v — 50. as in Fig. 5, c\ and C2 as in Eq. (101). 



fa = (1,0,0,-1), 
I = uJi(l,s h 0,ci), 

k = u k (l, s k cos 4>, s fc sin0, c k ). (114) 

Here we define sj^ = sin 6^ and = cos 6*^. 0i is the angle between the vectors / and Pi, 9 k is 
the angle between the vectors k and f3\ and <fi is the azimuthal angle of the gluon with momentum 
k relative to the plane defined by fli, $2 and I. The available phase space in polar angle for the 
radiated gluons is 6 k G (vr/2 - 5, tt/2 + 6) and B x G (0, vr/2 - 5) U (vr/2 + 5, vr). 

Using the diagrammatic rules for eikonal lines and vertices, as listed for example in [26], we 
can write down the expressions corresponding to each diagram separately. For example, diagram 
3 a) gives 

a) + (* « I) = [/*TX«AtJ] (-tfffflfff) W* + -0 ^7 

+ (A; <-> /). (115) 

V a/ 3^(k + Z, — /c, — /) = [(2& + l)^g a p + i} — ^)a9p~i — (2/ + ^/jflWy] is the momentum-dependent 
part of the three gluon vertex. Using the color identity / & c Tr (i a i&t c ) = iCfMcC 'a/2, and the 
approximation j3j ■ I ^> f3j ■ k for j = 1,2, which is valid due to the strong ordering of the final 
state gluon energies, we arrive at 

„) + (* „ o = ±<VA&C*rf^f L-jL-i + jj-lj-A . (116) 
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We proceed in a similar manner for the rest of the diagrams. The results are: 
b) + (k^l) ■■ 



1 n Kf n tPi-fo 2 1 



c) 
d) 

e) 

/) + (k <-> I) 



C F N c {C F -C A /2)g, 



k-l fc-kfo-V 
4 (A • P2) 2 1 



fc-lfo-lfc-kfo-k' 



117) 



The color factors in the last two equations of (117) are obtained from the identity Tr (t a t b t a t b ) = 
C f Mc{C f —Ca/2). Combining the terms proportional to the color factor C F J\fcC A , and including 
the complex conjugate diagrams, we find for the squared amplitude 



\M\ 2 = 2g i s C F N c C A !3 1 -f3 2 



+ 



Pi -fa 



• (US) 



Having determined the amplitude, we need to integrate |M| 2 over the phase space corresponding 
to the geometry given in Fig. 2. Specifically, we have to evaluate: 



1= — [dee-" 1 f ^ [ 
Af c J in (2tt) 3 2^ k 



d 3 / 



n (27r) 3 2a; fc J a (2tt) 3 2o; ( 



8{e-u k /y/s)8{e-fil,a))\M\ 2 , 



(119) 



where the weight function f(l,a) is given, as in Eqs. (4) and (11), by 

f(l,a) = 



^(l- Ci y- a st : 0,e(O, n/2-6) 
^(l + dY^st : 0,e(7r/2 + (J,7r), 



(120) 



with a < 1. 

Using the equalities: ft • (3 2 = 2, ft • I = u t (l - q), ft ■ I = cq(l + q), ft • A; = u k (l - c k ), 
ft ■ k = uj k {l + Ck) and k ■ I = uj k 0Ji{l — c k c\ — s k s\ cos0) in Eq. (118), performing the integration 
over 0, and changing the integration variable q — > — q in the angular region ^ G (tt/2 + 5, 7r), 
we easily arrive at the following three-dimensional integral: 



C F C A 
1 



a \ 2 1 /-sin 5 



7T / £ J — sin (5 
1 



f n5 dc fe ^ dq /^^e-^d-O 1 -^ 

J— sin (5 Jsin5 J Sxfs LU] 

1 



Cfc + Q 1 + C fe V 1 + Q ' I - C k J S k 1 + Ci 



1 1 



(121) 



We are interested in the (l/e)ln(l/e) behavior of /. This is obtained after performing the cq 



integral with the replacement e 



-^(l-Ci) 1 "^?/^ 



9(1 — — q) 1 a sf/y/s). Remainders do 
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not contain terms proportional to \ne. In this approximation, the q integration can be carried 
out, and we obtain the integral representation for the term containing ln(l/e): 



2 Cj?C a 



In 



ev ) 



15 dc k 



In 



s| — cos 2 8 / 
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1 + sin 5 



In 



1 + sin <T 



1 — sin 5 



(122) 



The potential non-global logarithm of e is replaced by ln(ev). The angular integral over Ck can 
be expressed in terms of dilogarithmic functions. The final expression for the term proportional 
to ln(ev)/e takes the form: 



C f Ca 



«,\ 2 1 
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In 
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(123) 



Equivalently, we can express our results in terms of the rapidity width of the region fl, Eq. (14), 
and we obtain 



/ = i-V-h.i 



7T 



ev 
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Ar] 



— + Ar] ( - In (2 sinh(A^)) ) + Li 2 
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The coefficient 



2 cos 



h(A^)J - Li2 (- 2 Bmh(A,/2) e^ 2 ) - Li 2 (e-^) 



C(Arj) = - - 



7T 



(124) 



(125) 



C F C A Hev) 

as a function of is shown in Fig. 8. Naturally, C is a monotonically increasing function of 
Ar/. For Ar/ — > 0, 

C ~ 0(Ar] InAr/) , (126) 

and the cross section vanishes, as expected. On the other hand, as the size of region Q increases, 
C rapidly saturates and reaches its limiting value [17] 



lim C = — 



(127) 



B Recoil 

In this appendix, we return to the justification of the technical step represented by Eq. (23). 
According to this approximation, we may compute the jet functions by identifying axes that 
depend only upon particles in the final states Nj c associated with those functions, rather than 
the full final state N. Intuitively, this is a reasonable estimate, given that the jet axis should be 
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Figure 8: C(Arj), as denned in (125), as a function of rapidity width Arj of the region Q. The 
dashed line is its limiting value, C(Ar] — > oo) = ir 2 /6. 

determined by a set of energetic, nearly collinear particles. When we make this replacement, how- 
ever, the contributions to the event shape from energetic particles near the jet axis may change. 
This change is neglected in going from the original factorization, Eq. (21), to the factorization 
in convolution form, Eq. (27), which is the starting point for the resummation techniques that 
we employ in this paper. The weight functions f N (Ni,a) in Eq. (21) are defined relative to the 
unit vector hi corresponding to a = 0, the thrust-like event shape. The factorization of Eq. (21) 
applies to any a < 2, but as indicated by the superscript, individual contributions to f N (Ni,a) 
on the right-hand side continue to depend on the full final state N, through the identification of 
the jet axis. 

To derive the factorization of Eq. (27) in a simple convolution form, we must be able to 
treat the thrust axis, hi, as a fixed vector for each of the states N s , Nj c . This is possible if 
we can neglect the effects of recoil from soft, wide-angle radiation on the direction of the axis. 
Specifically, we must be able to make the replacement 

fl(N Jc ,a)^f c (N Jc ,a), (128) 

where f c (Nj c , a) is the event shape variable for jet c, in which the axis h c is specified by state Nj c 
only. Of course, this replacement changes the value of the weight, e, f^(Nj c ,a) ^ f c (Nj ,a). 
As we now show, the error induced by this replacement is suppressed by a power of e so long 
as a < 1. In general, the error is nonnegligible for a > 1. The importance of recoil for jet 
broadening, at a = 1, was pointed out in [8]. We now discuss how the neglect of such radiation 
affects the jet axis (always determined from a = 0) and hence the value of the event shape for 
arbitrary a < 2. 

The jet axis is found by minimizing f(a = 0) in each state. The largest influence on the axis 
h c for jet c is, of course, the set of fast, collinear particles within the state Nj c associated with 
the jet function in Eq. (21). Soft, wide-angle radiation, however, does affect the precise direction 
of the axis. This is what we mean by 'recoil'. 
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Let us denote by uj s the energy of the soft wide-angle radiation that is neglected in the 
factorization (27). Neglecting this soft radiation in the determination of the jet axis will result 
in an axis hi(Nj c ), which differs from the axis fii(N) determined from the complete final state 
(N) by an angle A s 0: 

UMn)MNj c )) = A ^~^- ( 129 ) 

At the same time, the soft, wide-angle radiation also contributes to the total event shape 
f(N,a) ~ (l/Q)k a ± (k-y- a at the level of 

(130) 

because for such wide-angle radiation, we may take k~ ~ k S; ± ~ u s . In summary, the neglect of 
wide-angle soft radiation rotates the jet axis by an angle that is of the order of the contribution 
of the same soft radiation to the event shape. 

In the factorization (27), the contribution of each final-state particle is taken into account, 
just as in Eq. (21). The question we must answer is how the rotation of the jet axis affects these 
contributions, and hence the value of the event shape. 

For a wide-angle particle, the rotation of the jet axis by an angle of order A s in Eq. (129) 
leads to a negligible change in its contributions to the event shape, because its angle to the axis 
is a number of order unity, and the jet axis is rotated only by an angle of order e s . Contributions 
from soft radiation are therefore stable under the approximation (23). The only source of large 
corrections is then associated with energetic jet radiation, because these particles are nearly 
collinear to the jet axis. 

It is easy to see from the form of the shape function in terms of angles, Eq. (11), that for any 
value of parameter a, a particle of energy Ui at a small angle 9i to the jet axis fi\{N) contributes 
to the event shape at the level 

ei~|ft 2 -°. (131) 

The rotation of the jet axis by the angle A s due to neglect of soft radiation may be as large as, 
or larger than, 0j. Assuming the latter, we find a shift in the of order 



2-a 

Q V_SYV Q \ Q) Q 



Sii = e t (MAO) " Si (ni(N Je )) ~ £ (A s 0)^ a ~ £ I =f ] ~ % e s 2 ~ a . (132) 




The change in is thus suppressed by at least a factor e s l ~ a compared to e s , which is the 
contribution of the wide-angle soft radiation to the event shape. The contributions of nearly- 
collinear, energetic radiation to the event shape thus change significantly under the replacement 
(23), but so long as a < 1, these contributions are power-suppressed in the value of the event 
shape, both before and after the approximation that leads to a rotation of the axis. For this 
reason, when a < 1 (and only when a < 1), the value of the event shape is stable whether or 
not we include soft radiation in the determination of the jet axes, up to corrections that are 
suppressed by a power of the event shape. In this case, the transition from Eq. (21) to Eq. (27) 
is justified. 
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